Quark Orbital Angular Momentum in the MIT Bag Model 



Matthias Burkardt and Abdullah Jarrah 
Department of Physics, New Mexico State University, Las Cruces, NM 88003-0001, U.S.A. 

Using the MIT bag model, we study the contribution from the gluon vector potential due to the 
spectators to the orbital angular momentum of a quark in the bag model. For a s — 0(1), this spec- 
tator contribution to the quark orbital angular momentum in the gauge-covariant Ji decomposition 
is of the same order as the non gauge-covariant quark orbital angular momentum and its magnitude 
is larger for d than for it quarks and negative for both. 



PACS numbers: 



■ I. INTRODUCTION 

o : 

Since the famous EMC experiments revealed that only a small fraction of the nucleon spin is due to quark spins 
^ . [Hi there has been a great interest in 'solving the spin puzzle', i.e. in decomposing the nucleon spin into contributions 
O ' from quark/gluon spin and orbital degrees of freedom. In this effort, the Ji decomposition Q 

^1- 



<-£h appears to be very useful, as not only the quark spin contributions Aq but also the quark total angular momenta 
Jq = \Aq + L z q (and by subtracting the spin piece also the the quark orbital angular momenta L z ) entering this 
O i' decomposition can be accessed experimentally, through generalized parton distributions (GPDs). The terms in ([T]) 
O , are defined as expectation values of the corresponding terms in the angular momentum tensor 



M 0x V = ^ gt£* g + ^ g t (f x idy q + rx(SxBj 2 (2) 



in a nucleon state with zero momentum. Here iD = id — gA is the gauge-covariant derivative. The main advantages of 
this decomposition are that each term can be expressed as the expectation value of a manifestly gauge invariant local 
operator and that the quark total angular momentum J z — ^Aq+L z can be related to generalized parton distributions 
(GPDs) and is thus accessible in deeply virtual Compton scattering and deeply virtual meson production and can 
also be calculated in lattice gauge theory. 

Recent lattice calculations of GPDs surprised in several ways Q. First, the light quark orbital angular momentum 
7— ( , (OAM) is consistent with L u w —Ld, i.e. L u + Ld~0, which would imply that J g w | • 0.7 represents the largest piece 
in the nucleon spin decomposition. Secondly, L u s» —0.15 and L u s» +0.15 in these calculations, i.e. the opposite 
. i-H , signs from what one would expect from many quark models with relativistic effects, as we will also illustrate in the 
following section. While the inclusion of still-omitted disconnected diagrams may change the sum L u + Ld, it docs 
not affect the difference L u — Ld- In Ref. Q, it was pointed out that evolution from a quark model scale of few 
hundred MeV to the lattice scale of few GeV could possibly account for the difference. However, due to the presence 
of interactions through the vector potential in the gauge covariant derivative L z q does not have a parton interpretation, 
which complicates its physical interpretation. 

Jaffe and Manohar have proposed an alternative decomposition of the nucleon spin, which does have a partonic 
interpretation [j| 

^He a ?+e^ ag +^ (3) 

and whose terms are defined as matrix elements of the corresponding terms in the +12 component of the angular 
momentum tensor 



M' 



~ 12 = \ E 1+ ™+ +Y.1+( rfx ^) 2 1+ + £ + ~ ijr RF +i A j + 2TrF +j (r x id) * A 3 . (4) 



The first and third term in ([SI),© are the 'intrinsic' contributions (no factor of rx) to the nucleon's angular momentum 
J z = +| and have a physical interpretation as quark and gluon spin respectively, while the second and fourth term 



2 



can be identified with the quark/gluon OAM. Here q+ = 57 7+5 is the dynamical component of the quark field 
operators, and light-cone gauge A + = A + A z = is implied. The residual gauge invariance can be fixed by 
imposing anti-periodic boundary conditions A^(x^,oo~) = — Aj_(x^, — oo~) on the transverse components of the 
vector potential. 

Only the Aq are common to both decompositions. While for a nucleon at rest the difference in the Dirac structure 
between L z and C q plays no role 01 > the appearance of the gluon vector potential in the operator defining L z q implies 
that in general L z q 7= L z . Indeed, in Ref. |7j their difference was illustrated in the context of an electron in QED to 
order a. This lowest order correction to the electron orbital angular momentum from its own vector potential, can 
be directly translated into a QCD correction to the orbital angular momentum of a quark, resulting in [3] 

for a quark state with J z — +^ (and p = 0). However, in addition to this contribution to q^ (fx gA^j q from the A 

field 'caused' by its own current, the vector potential can also have been caused by the spectators, or may be due to 
intrinsic glue. In this note, we use the bag model to make an estimate for such 'spectator effects', by examining the 
contribution of the vector potential arising from spectator currents to L z — C z in perturbation theory to 0(pt s ). 



II. ORBITAL ANGULAR MOMENTUM IN THE MIT BAG MODEL 

In the bag model @, the Dirac wave function for a quark state with j z — s, is of the form 

^,=A^..^f;> ) ( 6) 

\ i]i(kr)r ■ aXs J 

where \s is a two-component Pauli spinor, k = E — 2, °^ 28 , with R being the bag radius, and N = J\f s is a normalization 
constant. For simplicity, we take the quarks to be massless. 

Evaluating the orbital angular momentum is straightforward in the absence of QCD corrections. For a state with 
s = +5 one finds from the lower component of the Dirac wave function 

C\ = |AA| 2 y J R drr 2 jl(kR) = 0.1735 (7) 

One can easily understand the sign of (JTJ) : since the quark spin s z = ± ^ and its orbital angular momentum must add 
up to j z — ^ for each wave function component, the ^-component of the quark orbital angular momentum can only 
be +1 or 0, i.e. < CZ < 1. For a quark state with j z — — ^, one obtains of course the opposite sign for C q . 
Using standard SU(6) wave functions, 

\p t) = — L {2 \u t u t d I) + 2 \u t d I u t) + 2 \d | u t u t) - \u t u | d f) - \u I u t d f) (8) 
v 18 



- |u t d t u \) - \u I d t u t) - \d t u t u X) - \d t u J, u f)} , 
this result gets multiplied by | for u quarks in a proton and by — ^ for d quarks in a proton, yielding 

C z u/p = ^-0.1735 = 0.2313 (9) 



C z d/p = ~. 0.1735 =-0.0578, 



with signs opposite to that of the lattice calculation [3J. A similar pattern for the signs as in Eq. © is observed in 
all quark models where the orbital angular momentum arises as a relativistic effect from the lower Dirac component. 
Similarly, phenomenological models for single-spin asymmetries have the same pattern of signs for quark angular 
momenta as the bag model (J9j) and thus also differ from the lattice QCD results Q. Although present lattice 
calculations still suffer from uncontrolled systematic errors due to the omission of operator insertions into disconnected 
quark loops, this does not affect the isovector combination — £-d/p> which is positive in the model calculations 
but negative in lattice QCD. Q 2 evolution has been proposed as a possible explanation to account for this apparent 
discrepancy 0, , but in order to accomplish quantitative consistency, perturbative evolution equations need to be 
used very deep in the nonperturbative regime. 
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Since quark models have, to lowest order in a s , no vector potential, it makes perhaps more sense to identify the 
quark OAM from these models with C z rather than with the GPD-based L z q . In Ref. [7| the difference A z = C z — L z 
was calculated to order a for a single electron in QED and the result then also applied to a single quark in QCD. 
However, in QCD quarks are never alone and the question arises regarding the effects from 'spectator currents' on 
the orbital angular momentum of each quark. 

In order to address this issue, we will in the following focus on estimating 0(a s ) corrections to the difference 



A q =L*-£ q = {ql[rxgA) q). (10) 

The vector potential in (|10j) is calculated from the spectator currents, which are obtained by taking matrix elements 
in the corresponding ground state bag model wave functions. The vector potential resulting from these static currents 
is obtained by solving 

V 2 A Q (r) = -f{f) = - £ frfl {r)a^ s , (r) (11) 

s' 

for each color component a and where the summation is over the spectators (here we pick the gauge V • A a = 0, but 
to 0(a s ) the result is actually gauge invariant as we will explain below). 
The contribution from a spectator with j z — s' to (q(rx A) z q) thus reads 

^' = ~H J d ^ d ^'^(^s{r)^) { j^^A^) (12) 

where the factor — | arises from the color part of the matrix element. Note that A s / is independent of the angular 
momentum j z = s of the 'active quark', since tp\ (r)ipi(f) — ip_ ± (f f )ip_i(r). However, it depends on the spin of the 

2 2 2 2 

spectator since the orientation of the vector potential entering (fT0|) depends on the latter. For example, for s' = +i, 
one finds 

4(f0(r X ct)'M?) = W?h{kr')h{kr')2 — ' + VV ' (13) 

and hence 

All = -\a s I R (14) 
+ 2 3 

In = \Af\*[ d*r[ MmW+JlW]* ^^^^ - (15) 

Jr<R Jr'<R I 5 r I r 

oo It 1 

The angular integration can be easily done using the expansion j^jL, = YliZo 2ITT Em=-l ^+r^ m (^)^/ m *(^') an d 
noting that only I = 1 contributes in (fT51) . yielding 

Ir = h( 4 ^) 2 |^| 4 [* 1 dr I dr'r' 2 r 3 [j$(kr) + jf(kr)] ^ (fcr')ji(fcr') = 0.1165 (16) 
y Jo Jo r > 

which does not depend on the bag radius. For s' — — h the result is identical with a negative sign. 

If the active quark has s alligned with that of the proton the two spectators must have opposite s' and their 
contribution to A cancels, i.e. A is nonzero only in those wave function components where the active quark has 
s = — |, in which case both spectators have s' = +5- As a result, A^ p is equal to twice A^_i times the probability 

to find that quark flavor with s = — | (which is | for q = u and | for q — d, and hence 

K/p = l A l± = -\<*sIr = -0.052a s (17) 



with 



A^ /p - = -\a s I R = -0.104a s 



which is the main result of this paper. 
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For simplicity, we evaluated the vector potential in the V • A — gauge. However, our perturbative result is gauge 
invariant (to 0{a s )) at least in the subclass of all gauges where all color components are treated (globally) SU(3)- 
symmetrically. I such gauges, matrix elements of operators of the type q 1 X a TqA a , where T is some Dirac matrix, and 
A a is calculated to 0(a s ), are proportional to the matrix elements of the corresponding abelian operators. Therefore it 
is sufficient to establish gauge invariance of q'r x Aq for abelian fields. The key observation is that the MIT bag-model 
wave functions contain no correlations between the positions of the quarks. Therefore, after eliminating the color in 
this calculation and introducing abelian currents, A z , factorizes into the density of the active quark ipl(f)ip s (r) times 
(r x A) z = rA^. Writing the volume integral J d 3 r in cylindrical coordinates, one can isolate the only ^-dependent 

term r d<pA,p — § dr ■ A(f) as a closed loop integral with fixed r and z. The closed loop integral is gauge invariant 
(its numerical value represents the color-magnetic flux through a circle with radius r) and so is the volume integral 
in which it enters. 

In addition to ([T]) and ([31), other angular momentum decompositions have been proposed in the literature (see for 
example Refs. [1, Since we used Coulomb gauge, A = A p h ys , where A^hys is the transverse part of the vector 
potential as defined in Refs. [H, 0]. In the decomposition proposed in Ref. Q, the quark orbital angular momentum 
involves only the longitudinal part of the vector potential, which is zero in our calculation, and which, in the bag 
model, as well as many other quark models, vanishes also in other gauges to 0(a s ) as discussed above. Therefore, 
within our approximations, the quark orbital angular momentum as introduced in Ref. Q would not be affected by 
spectator effects and would be given by Eq. ©. In contradistinction, the quark orbital angular momentum in Ref. 
Q is identical to that of [2| and therefore would receive the same contribution from the vector potential as L z . 



III. DISCUSSION 



We found that the effect from the vector potential caused by spectators on A* , p = L z — £ z is negative for both 

q = u,d. This is consistent with the sign that one may have intuitively expected: since A is obtained using the 

Biot-Savard equation A(f) — J d 3 r' \ \ , one would expect that A and j have the same orientation. Hence r x A 

should have the same sign as the magnetization density rxj. In the sum over all spectators, quarks with j z = +i 
dominate and thus the magnetization density of the spectators has the same sign as the nucleon spin. Taking into 
consideration the negative color factor which reflects the fact that the interaction between the quarks is attractive, 
and the fact that only spectators with s' contribute, one arrives at the above negative sign. 

The difference A* , — A^ p — 0.052a s is positive. This implies that the vector potential from the spectators adds 
a negative contribution to L z u — C z d when compared to — L z d . Since — L z d calculated in lattice QCD is already 
negative, subtracting our result for A* , — ^d/p wou ld make it even more negative, and therefore this does not help 
in understanding the puzzling lattice results. 

The numerical values of A* , p are of the same order of magnitude as the corresponding difference for a free quark [3] ■ 
Whether our result is overall numerically significant depends on the numerical value of a s — a matter of debate in the 
context of the bag model. For example, in order to generate the experimentally measured N — A mass splittings from 
perturbative color-hyperfine interactions in the bag model, a s = 0(1) is required — in which case the A*, p would be 
large enough that L d / p would become positive, but still not large enough to render L u / p negative. Of course, such a 
large value of a s is inconsistent with the use of perturbation theory, and may only reflect the need for correlations 
between quark wave functions, which are not present in the MIT bag model. However, nonperturbative correlations 
that enhance the hyperfine splitting, may also enhance A z . Since both A^ and A d are negative, their difference is 
small which seems to indicate that spectator contributions to the isovector combination L u — L d are small, but may 
not be small in the case of the isoscalar combination L u + L d . 
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